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I.  INTRODUCTION 


To  properly  describe  chemical  reactions 
in  non-permeable  porous  media,  the  "pore  tree" 
was  introduced  by  Simons  and  Finson  (1979) 
and  Simons  (1982).  This  pore  structure  was 
developed  via  analogy  to  the  kinetic  theory  of 
gases;  the  pore  length  is  analogous  to  the  mean 
free  path  and  under  the  assumption  that  the 
pore  aspect  ratio  (pore  length  to  radius)  is  a 
constant,  a  pore  size  distribution  was  obtained 
that  has  been  confirmed  for  coal,  coal  char, 
sorbents,  catalysts  and  even  for  kidney  stones 
from  both  men  (low  porosity  oxalate)  and 
women  (high  porosity  phosphate).  The  pore 
tree  was  statistically  derived  from  the  pore  size 
distribution  and  allows  the  orderly  migration  of 
a  reactant  gas  from  the  large  pores  to  the  small 
pores  (see  Fig.  1).  A  detailed  description  of 
the  pore  tree  and  the  coupled  transport  and 
chemistry  is  given  by  Simons  (1982,  1983a). 

The  spatially  dependent  transport/reaction  equations  are  solved  for  a  single  pore  tree  and  then 
the  total  contribution  of  all  trees  (of  all  sizes)  in  the  system  is  obtained  by  summing  the 
contribution  of  each  tree  that  reaches  the  exterior  of  the  system.  This  is  distinct  from  the 
"bulk"  transport  equations  that  integrate  over  all  pores  at  a  fixed  point  in  space  before  they 
integrate  spatially.  The  "bulk"  transport  approach  inherently  assumes  instantaneous  mixing 
between  large  and  small  pores  at  a  single  point  in  space  (Random  Pore  Models)  while  the 
pore  tree  relaxes  this  assumption  with  a  structure  that  is  developed  statistically  and  is 
relatively  easy  to  work  with  mathematically. 

In  all  previous  applications  (sorbents,  catalysts  and  coal  combustion)  the  pore  tree 
represents  an  isolated  sub-structure,  allowing  diffusion  into  and  out  of  the  porous  media 
without  permitting  convection  through  the  media.  Under  the  current  ARO  program,  Simons 
Research  Associates  is  adapting  its  pore  tree  model  to  describe  the  permeable  pore  structure 
and  convective  transport  characteristic  of  the  subsurface  flow  of  water  in  soil,  the  dispersion 
of  contaminants  and  the  in-situ  remediation  of  contaminated  sites.  The  random  nature  of  the 
pore  structure,  which  formed  the  basis  of  the  statistical  derivation  of  the  pore  tree,  is  applied 
to  the  underground  structure  of  porous  soil,  sand  and  rocks.  The  statistical  description  of  the 
pore  tree  is  extended  to  describe  the  concept  of  permeability,  thus  allowing  convective  flow 
as  a  transport  mechanism  in  addition  to  the  diffusive  transport  already  considered.  Physically, 
this  amounts  to  statistically  determining  the  "branches"  that  are  common  to  several  trees  to 
allow  percolation  through  the  large  scale  (mobile)  structure  in  addition  to  diffusion  through 
the  smaller  scale  (immobile)  structure.  While  permeability  is  dominated  by  the  largest  pores, 
it  is  imjKwtant  to  determine  the  level  of  convection  that  is  occurring  at  the  intermediate  scales 
in  order  to  accurately  describe  transport.  The  resulting  pore  structure  will  provide  an  analytic 
description  of  an  underground  network  upon  which  transport  and  coupled  chemical  reactions 
may  be  accurately  superimposed. 


/yA  fp  *  Appropriate 
Mean  Radius  at  x 

k 

Particle  Surface 


Figure  1.  The  Pore  Tree 
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n.  ISOLATED  PORE  TREE;  THE  STRUCTURE 

Following  the  pore  structure  theory  of  Simons  and  Finson  (1979)  and  Simons  (1982), 
consider  a  spherical  porous  particle  of  radius  a,  containing  pores  of  length  Ip  and  radius  rp. 
The  pore  dimensions  range  from  a  microscale  of  the  order  of  Angstroms  to  a  macroscale 
which  is  a  significant  fraction  of  the  particle  radius.  The  radius  of  the  largest  pore  is  denoted 
by  r„„„  and  is  given  by 


where  6  is  the  total  porosity  of  the  particle  and  K„  is  a  constant  of  integration, 
approximately  equal  to  five,  which  relates  the  pore  length  to  its  radius 


The  radius  of  the  smallest  pore  is  denoted  by  r,„i„  and  is  given  by 


'•mm=2e/Pp/p 

where  p,  is  the  density  of  the  solid  matrix,  Sp  is  the  specific  internal  surface  area  (several 
hundred  m^/g),  and 


The  particle  contains  a  continuous  distribution  of  pore  sizes  from  r,„i„  to  .  The 
number  of  pores  within  an  arbitrary  plane  of  cross-sectional  area  A  and  with  radius  between 
rp  and  rp  +drp  is  denoted  by  g(r^dr^.  The  pore  distribution  function  g(rpis  given  by 

g(/-p.e/2;lPr; 

where  indicates  an  average  over  all  inclination  angles  between  the  axis  of  the  pore  and 
the  normal  to  the  plane.  Due  to  the  random  orientation  of  the  pores,  the  intersection  of  a 
circular  cylinder  with  a  plane  is  an  ellipse  of  average  area  2TZTp.  Hence,  the  porosity  is  the 
27trp^  moment  of  g(r^  and  the  internal  surface  area  is  the  47irp  moment  of  g(r^ .  The 
expression  for  g(r^)  was  derived  (Simons  and  Finson,  1979)  from  statistical  arguments  and  has 
been  validated  through  extensive  comparison  of  the  predicted  volume  and  surface  area 
distributions  with  mercury  intrusion  data  (Stacy  and  Walker,  1972).  This  has  been 
accomplished  for  coal,  char  derived  from  that  coal  (Kothandaraman  et.al.,  1984),  sorbents, 
catalysts  and  even  kidney  stones  from  both  men  (oxalate)  and  women  (phosphate). 
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A  characteristic  feature  of  the  1/rp^  distribution  depicts  that  the  pore  volume  between 
r  .  and  tp  increases  linearly  with  the  natural  log  of  rp.  It  is  the  functional  form  of  this 
relationship, 


Pore  Volume 


r 


Inr. 


*Tntn 


(6) 


that  depicts  the  inverse  cubic  dependence  of  g(r^  on  rp.  A  linear  display  of  mercury  intrusion 
volume  vs.  In(rp)  always  infers  a  l/Vp^  distribution. 

The  number  of  f)ores  within  the  bulk  volume  V  whose  pore  radius  is  between  and 
rp+dtp  may  be  defined  by  Vfir^r^.  The  pore  volume  is  expressed  as  the  Jtrp\  moment  of 
/(rp  and  the  internal  surface  area  is  the  27t:rplp  moment  of  /(rp.  The  pore  size  distribution 
fiinctions  (  /(rp  and  g(r p  )  are  clearly  not  independent.  The  definitions  of  porosity  and 
internal  surface  area  infer  that  /(rp  is  related  to  g(r^  by 

Equation  (7)  simply  states  that  the  probable  number  of  pores  intersecting  an  arbitrary  plane 
increases  with  the  length  of  the  pore  and  with  the  density  of  pores. 

The  length  of  a  pore  is  determined  by  an  arbitrary  intersection  with  another  pore  and 
is  expressed  (Simons  and  Finson,  1979)  as  a  collision  integral  over  the  pore  distribution 
functions.  The  analysis  suggests  that  Ip  ,  ^(rp  and  /(rp  are  proportional  to  rp  ,  I/t^^  and 
respectively.  The  constants  of  proportionality  are  obtained  from  integral  constraints,  i.e.,  the 
total  porosity  and  internal  surface  area  contained  in  the  pore  structure.  The  expression  for 

Arp  is  given  by 


where  the  constants  were  defined  above. 

The  pore  volume  distribution  corresponding  to  these  distribution  functions  is  similar  to 
that  utilized  in  the  random  pore  model  (Gavalas,  1980  &  1981).  However,  the  pore  tree 
model  and  the  random  pore  model  differ  dramatically  in  their  choice  of  the  pore  aspoct  ratio 
(length  to  diameter)  and  its  implications  with  respect  to  pore  branching.  The  random  pore 
model  allows  a  single  pore  to  connect  two  larger  pores.  This  picture  lends  itself  to  the 
idealization  of  instantaneous  mixing  between  the  pores  and  requires  that  the  pore  aspoct  ratio 
be  of  the  order  of  one  hundred.  The  pore  tree  theory  uses  data  for  r,„^  to  imply  (via  K„)  that 
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all  pores  possess  an  aspect  ratio  of  the  order  of  ten.  Hence,  small  pores  may  connect  to  larger 
pores  only  on  one  end  and  all  pores  must  branch  from  successively  larger  pores  like  a  tree  or 
river  system. 

Each  pore  that  reaches  the  exterior  surface  of  the  particle  is  depicted  as  the  trunk  of  a 
tree.  The  size  distribution  of  tree  trunks  on  the  exterior  surface  of  the  particle  is  denoted  by 

where  g{r^  is  functionally  identical  to  g{r^.  Each  trunk  of  radius  r,  is 
associated  with  a  specific  tree-like  structure.  Let  N,  be  defined  as  the  branch  distribution 
function  where  Njdrp  is  the  number  of  pores  of  radius  rp  (within  size  range  drp)  in  a  tree 
whose  trunk  radius  is  ri .  The  total  number  of  pores  of  radius  rp  in  a  sphere  of  radius  a 

may  be  expressed  as  4/3Jta^y(rp</r^  or,  as  the  sum  of  all  pores  of  radius  rp  contained  within 
every  tree  in  the  porous  sample,  plus  all  pores  of  radius  rp  that  are  themselves  the  trunk  of  a 
tree.  Hence, 


^Tia^Ar^=  f  N^(r,)4na^dr^  +  4na^g(r^ 


(9) 


where  g(r^)  is  the  number  of  tree  trunks  per  unit  external  area  of  the  porous  sample  and  only 
those  trees  whose  trunk  radius  is  greater  than  rp  may  contain  a  pore  of  radius  rp .  Using  the 
previously  derived  expressions  for  r^,  g^r^and  A''p)y  E^-(9)  is  identically  satisfied  by 


(10) 


The  branch  distribution  function  completely  characterizes  the  pore  tree.  The  internal 
surface  area  and  pore  volume  associated  with  each  pore  tree  are  denoted  by  St(rJ  and  Vj(rJ, 
respectively,  and  are  expressed  as  the  sum  of  the  contributions  from  the  trunk  and  that  from 
the  branches. 


S,(0=2nr/  +  f 


(11) 


(12) 


'  ada 


Using  Eq.(lO)  for  N„  and  V^rJ  become 


5 


(13) 


\  min  y 


n-6) 


\  V^minjy 


(14) 


where  the  (1-0)  term  in  St  has  been  included  to  account  for  pore  combination  (Simons, 
1979a). 


The  surface  area  associated  with  the  pore  tree  may  be  several  orders  of  magnitude 
greater  than  the  surface  area  of  the  trunk.  However,  the  volume  of  the  pore  tree  may,  at 
most,  be  one  order  of  magnitude  greater  than  that  of  the  trunk.  It  should  also  be  noted  that 
the  above  expressions  for  Sj  and  V,  reduce  to  those  appropriate  to  a  single  cylindrical  pore  in 
the  limit  of  rj  (the  leaf  of  the  tree).  Furthermore,  the  integrals  of  St(rJ  and  Vt(rJ  over 
aU  g(r)  recover  the  total  internal  surface  area  and  pore  volume  of  the  porous  sample. 

Each  trunk  of  radius  r,  is  associated  with  a  specific  tree-like  structure  with  continuous 
branching  to  ever  decreasing  pore  radii.  The  radius  and  number  of  pores  is  a  unique  function 
of  the  distance  x  into  the  tree.  The  coordinate  x  is  skewed  in  that  it  follows  a  tortuous  path 
through  the  branches  of  the  tree.  Let  n(x)  represent  the  number  of  pores  of  radius  rp  at 
location  x  in  a  tree  of  trunk  radius  r,.  An  analysis  (Simons,  1982)  of  this  pore  tree  has 
demonstrated  that 


nix)=rf/r^(x) 

and  the  coordinate  x  is  related  to  rp  by 

drjdx=-rjl^ 

The  continuous  branching  model  has  been  used  to  successfully  describe  char  oxidation 
(Lewis  and  Simons,  1979;  Simons,  1982  &  1983a),  coal  pyrolysis  (Simons,  1983b  &  1984) 
and  the  catalytic  cracking  of  benzene  by  porous  iron  oxides  (Simons  et  al.,  1986).  It  was  also 
used  to  successfully  describe  sulfur  sorption  (SO,  and  H2S)  by  porous  calcine  (CaO)  in  the 
limit  of  zero  utilization  (Simons  and  Rawlins,  1980;  Simons  et  al.,  1984)  and  was  later 
extended  to  include  CaS04  and  CaS  deposits  (Simons  and  Carman,  1986;  Simons  et  al.,  1987; 
Simons,  1988;  Simons  et  al.,  1988).  The  subsequent  determination  of  the  controlling  physical 
parameters  led  to  a  new  concept  for  the  optimization  of  the  sulfur  sorption  process  (Simons, 
1991;  Simons  et  al.,  1992)  through  spray  drying  of  water  soluble  organic  calcium  solutions  to 
control  the  sorbent  pore  structure. 


(15) 


(16) 
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III.  INTERCONNECTIVITY  AND  PERMEABILITY 

The  first  step  in  determining  the  size  distribution  of  the  interconnected  pores  and  the 

distribution  of  the  permeability  is  to  determine  the  distribution  function  Gf(r^,r^dr^  which 
represents  the  number  of  pores  of  radius  (within  size  range  drp)  per  unit  cross  section  of  an 
arbitrary  plane  and  also  contained  within  a  tree  whose  trunk  radius  is  r,.  Consider  an  infinite 
homogeneous  isotropic  porous  .media  and  isolate  a  spherical  volume  of  that  media  denoted  by 
the  radius  a.  Such  a  volume  is  illustrated  in  Fig.  2.  The  total  number  of  pores  of  radius  rp 
(within  size  range  drp)  intersecting  plane  AA  of  area  Tia^  has  previously  been  defined  by 

g{r^%a^drj^  .  The  pores  in  plane  A  A  in  this  size  range  may  also  be  determined  by 

integrating  Gj(r^^r^%a^drp  over  all  trees  whose  trunk  intersects  the  exterior  surface  of  the 
porous  sample.  Hence  it  follows  that 


g(r^ita^dr^=  J  [G,(r^yr^na^dr^g(r;)4na^dr^ 


(17) 


where  only  those  trees  whose  trunk  radius  is  greater  than  rp  may  contain  a  pore  of  radius  rp . 


A  solution  to  Eq.  (17)  for  Gr,(r,,rp  will  not  necessarily  be  unique.  Physical  arguments 

will  help  determine  G,(rj,rp  and  help  ensure  that  it  is  the  particular  solution  we  seek.  Since 
N,  represents  the  number  of  pores  of  size  rp  in  the  tree  and  the  probability  of  a  pore 

intersecting  a  plane  is  proportional  to  its  length,  it  follows  that  Gf(r^,r^  should  be 
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proportional  to  the  product  of  Nj  and  XJ\  ,  i.e.,  proportional  to  Eq.  (17)  is  identically 

satisfied  by  a  function  which  differs  from  t^lx^  by  In(rp). 


(18) 


Note  that  the  ln(rp/r^)  introduces  an  integrable  singularity  at  rp=r^.  No  physical 
significance  to  this  singularity  has  yet  been  identified. 


Further  confirmation  of  the  form  for  may  be  made  through  the  definitions  of 

porosity  and  internal  surface  area.  If  the  porosity  and  internal  surface  area  in  plane  AA  due  to 
a  single  tree  of  radius  q  are  defined  by 


'  -  r 

0,  =  /  0,(r„r,)2xr^dr^  » 


(19) 


and 


s,  =  f  —G,(r,,r^)4%r^dr^  « 

r  ■  Ps 


(20) 


respectively,  then  the  total  porosity  and  internal  surface  area  must  be  represented  by  their 
respective  integrals  of  0t  and  St  over  all  trees.  Hence  6 

^nmx 

0  =  f  0,g(rt)4na^dr^ 

and  Sp 

^amx 

=  f  s,g(r)4na^dr^ 

^rafa 

follow  respectively. 

It  is  stated  without  further  proof  that  Eqs.  (19)  to  (22)  are  compatible  in  both  their 

exact  and  their  approximate  forms.  The  distribution  function  Gf(r^,r^)  as  given  by  Eq.  ( 1 8) 
is  exact  and  will  be  used  to  determine  the  interconnectivity  of  the  pore  structure. 
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The  probability  of  trees  sharing  common  branches,  i.e.,  the  interconnectivity  of  the 

pore  structure  is  described  in  Fig.  3.  We  seek  the  distribution  function  which 

represents  the  number  of  pores  (within  size  range  drp  about  rp)  per  unit  area  of  plane  AA  that 
are  connected  to  both  sides  of  the  pore  structure  through  pores  at  least  as  large  as  rp.  A^  is 
defined  as  the  area  within  plane  AA  that  is  open  to  one  side  of  the  porous  media  through  all 
trees  of  size  r,'  (through  all  pores  of  size  rp’  that  are  at  least  as  large  as  rp).  Subsequently, 

^^®  humber  of  pores  of  size  rp  (within  size  range  drp)  per  unit  area 
of  plane  AA  that  are  contained  in  a  tree  of  size  range  dr,  about  r,  and  are  also  connected  to 
the  opposite  side  of  the  porous  media  through  all  trees  denoted  by  r,'.  It  follows  that  the 
distribution  function  for  interconnected  pores  in  plane  AA  may  be  obtained  by  integrating 

\0,(r„r^dr^  over  all  trees  (rj  that  are  large  enough  to  contain  a  pore  of  size  rp  .  Hence, 


Jir^%a^drp=  f  ^e^,(r,,r^dr^g(r,)2na^dr^ 

'> 

From  the  above  definition  of  A^  ,  A^  may  be  expressed  as 

f[f  2itrpKa^G,ir,,r^dr^g(r)2ita^dr^ 


where  the  primes  on  the  variables  of  integration  have  been  omitted.  Evaluating  Eq.  (24) 


Figure  3.  Interconnectivity  of  a  Porous  Media 


(23) 


(24) 
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yields 


TCfl^ein(r^/rp 

2p 


(25) 


from  which  Eq.  (23)  yields  the  common  branch  distribution  function. 


ehl(r 


4P 


—  s(V 


(26) 


It  has  been  deduced  that  the  total  number  of  common  branches  of  size  rp  in  an 
arbitrary  plane  scales  approximately  with  the  total  number  of  pores  of  that  size  in  that  plane. 
Eq.  (26)  suggests  that  this  proportionality  is  logarithmic  in  pore  size  and  that  approximately 
one  percent  of  all  pores  of  all  sizes  are  interconnected  through  larger  pores.  Hence, 
convection  will  be  responsible  for  transport  in  approximately  one  percent  of  the  small  pores. 
While  permeability  is  dominated  by  the  largest  pores,  it  is  important  to  determine  the  level  of 
convection  that  is  occurring  in  smaller  pores  in  order  to  accurately  describe  the  fine  scale 
transport  necessary  to  assess  chemical  reactions.  The  level  of  interconnectivity  deduced  above 
suggests  that  a  very  wide  range  of  pore  sizes  separates  the  "mobile"  and  "immobile"  regions 
of  the  pore  structure  and  that  a  complicated  mixture  of  convective  and  diffusive  transport 
persists  in  this  pore  size  range.  Determining  permeability  as  a  function  of  pore  size  is  the  first 
step  in  describing  the  transport  in  this  pore  size  range. 

Permeability  across  plane  AA  will  possess  contributions  from  three  source:  first,  the 
large  tree  in  Fig.  2  may  provide  a  direct  link  across  the  porous  sample;  second,  the 
interconnectivity  derived  above  will  provide  that  contribution  from  pores  of  size  rp  that  are 
interconnected  through  larger  pores;  and  third,  a  pore  of  size  rp  may  be  interconnected 
through  smaller  pores.  Imagine,  for  example,  the  large  tree  in  Fig.  2  whose  fine  branches 
intersect  another  tree  rather  than  to  directly  exit  the  porous  sample.  This  connectivity  would 
translate  into  a  slower  velocity  in  the  larger  pore  (rp)  crossing  plane  A  A.  Each  of  these 
sources  will  be  evaluated  in  subsequent  technical  reports. 
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IV.  SUMMARY 


This  report  is  the  second  technical  report  on  the  project  entitled  "Development  of  Pore 
Structure  Models  for  Water  and  Contaminant  Transport  in  Partially  Frozen  Soils,"  submitted 
to  ARO  for  the  period  13  October  1994  to  31  October  1994.  The  pore  tree  model  is  being 
extended  to  describe  the  permeable  pore  structure  which  characterizes  the  subsurface  flow  of 
water  in  soil,  the  dispersion  of  contaminants  and  the  in-situ  remediation  of  contaminated  sites. 
Permeability  requires  a  statistical  determination  the  "branches"  that  are  common  to  several 
trees  to  allow  percolation  through  the  large  scale  (mobile)  structure  in  addition  to  diffusion 
through  the  smaller  scale  (immobile)  structure.  While  permeability  is  dominated  by  the  largest 
pores,  it  is  important  to  determine  the  level  of  convection  that  is  occurring  at  the  intermediate 
scales  in  order  to  accurately  describe  transport.  Approximately  one  percent  of  all  pores  are 
interconnected  through  larger  pores.  Hence,  convection  will  be  responsible  for  transport  in 
some  of  the  small  pores  without  influencing  the  over  all  permeability.  This  suggests  a  very 
wide  range  of  pore  sizes  separating  the  mobile  and  immobile  regions  and  a  complicated 
mixture  of  convective  and  diffusive  transport  in  this  pore  size  range.  Subsequent  analysis  will 
utilize  the  interconnectivity  of  the  pores  to  determine  the  distribution  of  the  permeability  with 
pore  size. 
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